Theory for the excitation spectrum of High-T c superconductors 
quasiparticle dispersion and shadows of the Fermi surface 
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Using a new method for the solution of the FLEX-equations, which allows the determination 
of the self energy Ek(w) of the 2D Hubbard model on the real frequency axis, we calculate the 
doping dependence of the quasi-particle excitations of High-T c superconductors. We obtain new 
results for the shadows of the Fermi surface, their dependence on the deformation of the quasi 
particle dispersion, an anomalous w-dependence of ImEk(<^) and a related violation of the Luttinger 
theorem. This sheds new light on the influence of short range magnetic order on the low energy 
excitations and its significance for photoemission experiments. 
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Despite an enormous progress jjj, the electronic exci- 
tation spectrum of the strongly correlated High-T c su- 
perconductors is still far from being understood. Re- 
cent angular resolved photoemission measurements 
for different doping concentrations demonstrate that pro- 
nounced deformations of the quasiparticle dispersion oc- 
cur. The opening of a spin density gap and the variation 
in spectral weight as function of k || reflect the strong in- 
fluence of the antiferromagnetic correlations on the low 
energy excitations. In particular, the interpretation of 
the shadows of the Fermi surface (FS) observed by Aebi 
et al. [Q in terms of antiferromagnetic correlations is un- 
der current debate [Q. Therefore, the determination of 
the elementary excitations with high energy and momen- 
tum resolution is of extreme importance. 

In this Letter, we calculate the quasiparticle excitation 
spectrum of the one band Hubbard Hamiltonian using 
a new numerical method for the self consistent summa- 
tion of all bubble and ladder diagrams g (fluctuation 
exchange approximation) on the real frequency axis, and 
compare our results with recent photoemission experi- 
ments. Due to the high degree of numerical stability 
of this method, we present interesting new results for a 
strong deformation of the quasiparticle dispersion at the 
wave vector k = (tt, 0) upon doping, the occurrence of 
satellite peaks and shadows of the Fermi surface in the 
paramagnetic state and an unusual momentum and fre- 
quency dependence of the electronic self energy. All this 
sheds new light on the occurrence of FS-shadows w ithout 
long range antiferromagnetic order. 

The fluctuation exchange (FLEX) approxima- 
tion &jllj is conserving in the sense of Kadanoff and 
Baym |12fl . It is a theoretical approach complementary 
to the exact diagonalization studies |l3] or quantum 
Monte Carlo simulations |l4|,[l5|], which gave deep in- 
sight into the origin of the quasiparticles of strongly 
correlated systems, but are limited to rather small sys- 
tems. The self energy in the FLEX approximation of 
the paramagnetic phase is given by the momentum and 



Matsubara-frequency sum = ^2 k , V(k — k')G(k'), 
where k = (k,iu) n ). The effective i nteraction V{q) re- 
sults from the summation of bubble and ladder diagrams 
and can be expressed in terms of the particle-hole bub- 
ble X (q) = -wE fc G(fc + q)G(k). § G(k) is the fully 
renormalized Greens function, uj n a fermionic Matsubara 
frequency, k the crystal momentum, N the number of 
lattice sites, and T the temperature. 

Because of the straightforward applicability of the 
numerically very effective fast Fourier transformation 
(FFT) for the evaluation of the above convolutions, most 
of the solutions of the FLEX equations are performed on 
the imaginary frequency axis. Unfortunately, from the 
imaginary axis, reliable information about the real fre- 
quency excitations can only be obtained using rather ill 
defined numerical procedures. Based on the contour inte- 
gral technique, we present in the following a numerically 
stable calculation of the electronic self energy on the real 
axis using the FFT. Note that recently Dahm et al. |ll]] 
also solved the FLEX equations on the real axis. How- 
ever, they performed the direct calculation of the rather 
time consuming frequency integrals. 

In the following, we first apply our approach to the 
one band Hubbard model in the paramagnetic state. The 
bare dispersion is given by e£ = — 2t (cos(k x ) +cos(k y )) — 
/i, with nearest neighbor hopping element t Q — 0.25 eV, 
and chemical potential fj,. The local Coulomb repulsion 
is given by U = 4t . The extension to multi band Hamil- 
tonians and to superconducting or magnetic phases is 
straightforward. |ll| We solve the set of equations for 
complex frequencies z = u) + i"f with small but finite 
imaginary part 7 < ttT. Then, the analytical continua- 
tion to the real axis 7 — > + is numerically well defined. 

The above momentum summations become simple 
products in Wannier space. Furthermore, the frequency 
summations are evaluated using the contour integral 
technique jLffl , where we place the four horizontal lines of 
the contour by a finite amount 7 away from the poles of 
the Greens function. Considering the particle hole bub- 
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ble Xi( z ) & t lattice site i and for frequencies z = uo + £7, 
we perform the Kramers Kronig transformation for the 
Greens function and express the resulting energy de- 
nominators using Laplace transformation [ jl8| . Now, the 
occurring frequency integrals decouple and we obtain 
Xi(z) = f™dtxi(t)e izt with 

Xi (t) = -i{2nfe-^ ( gi (t){AU(t) + A-ii-t^ 1 } 

- g *_.(t)[Mt)+At(-t)e 2 ^}) . (1) 

Here, A\{t) is defined by 

/OO J 
^r(e + i 7 )Gt(e + i 7 )e-^, (2) 
-00 27T 

and Qi(t) is the Fourier transform of the spectral density 
— ilmGi(e + £7). Note that we are using the complex 
Fermi function f(z), because we solve the FLEX equa- 
tions consistently for finite 7. This is essential for the fine 
structure at small energies. Thus Xk(z) and finally Vk(z) 
can be evaluated by a D + 1 dimensional FFT of Xi(i). 
The procedure for the frequency summation of the self 
energy is similar. Within the Wannier space, it follows 
for the Laplace transform of the self energy: 

Ei(t) = -i2irTVi(0) g,(t) + i^rfe^ (v l (t)[Ai(*)+ 

At{-t)e^] - ei (t)[Bt(t) +B i (-t)e 2 '< t }) . (3) 

Vi(t) is the Fourier transform of ——ImVi(z) and Bi{t) 
can be obtained from Eq. || by replacing f*{z)G\{z) by 
n* (z)V-* {z) . n(z) is the complex Bose function, and the 
first term results from the coincidence of the fermionic 
Matsubara frequencies uj n = oj n i , which cannot be in- 
cluded in the standard contour. Thus, the self energy 
Ek(w + 17) results via D + 1 dimensional FFT from 
Si(t) and the limit 7 — > + can be performed. 

The advantage of our method compared to the first 
real frequency calculation recently published by Dahm 
et al. is the applicability of the FFT and the con- 
sistent solution for a small but finite imaginary part 7, 
which lead to an accuracy of the convergence of the spec- 
tral density per momentum and frequency point of 10~ 5 . 
Only this accuracy enables us to obtain the new physical 
results presented here [ p0[ . 

In the inset of Fig. [I], we show results for the local 
(momentum averaged) density of states (DOS) for vari- 
ous doping values. We find for larger doping a rigid band 
like behavior, but for smaller doping a pseudogap occurs 
as a precursor of the Mott-Hubbard splitting. From the 
asymmetry with respect to the pseudogap, a transfer of 
spectral weight from high energy to low energy scales H] 
can be seen. 
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FIG. 1. Fermi surface in the quarter of the Brioulline zone 
with k x , k y > (solid line) in comparison with the Fermi sur- 
face of the uncorrelated system for different doping concen- 
trations. Note the violation of the Luttinger theorem and the 
occurrence of FS-shadows for smaller doping (n — 0.44). The 
inset shows the local density of states for various doping val- 
ues. Note, the occurrence of a pseudogap near the half filled 
case. 

In Fig. [I], results are shown for the FS for different 
band fillings in comparison with that for U — 0. The FS 
is obtained from the k-points, where e£ + ReEk(0) = 0. 
The changes of the FS are rather small, but a tendency to 
a stronger nesting topology is clearly visible. Such a de- 
formation of the FS is related to a dramatic deformation 
of the quasi particle dispersion. For larger doping, the 
Luttinger theorem |22j is fulfilled. However, for smaller 
doping the volume of the FS decreases compared to the 
U = case. This violation of the Luttinger theorem re- 
sults from a transfer of particles to additional shadow 
states (schematically shown in Fig. [l]), as will be dis- 
cussed below. 

In Fig. ||, we show our results for the quasi particle dis- 
persion obtained from the maxima of the spectral density. 
Since we expect the FLEX to be a good approximation 
for the low energy excitations, we focus on the dispersion 
near the Fermi energy and in particular in the neighbor- 
hood of the X-point. For larger doping (dashed line), 
the saddle point, responsible for the van Hove singular- 
ity, is visible. However, for smaller doping (solid and 
open-crossed squares) pronounced deformations of the 
dispersion occur. The dispersion of the states with high 
spectral weight (solid squares) is flattened and suddenly 
repelled from the Fermi energy. This might be a precur- 
sor of the large spin density gap at the X-point, found in 
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recent photoemission experiments by Wells et al. |p| for 
undoped cuprates. Furthermore, new quasiparticle states 
with low spectral weight (open-crossed squares) occur. 
These states, which are for larger doping only visible far 
away from the Fermi level (not shown), remain stable 
for low excitation energies and build up the shadows of 
the Fermi surface 0, which one usually would expect 
only for systems with long range antiferromagnetic or- 
der. In order to check this interpretation, we performed 
for the first time the spin resolved FLEX approximation 
within a given antiferromagnetic background (solid line) . 
We fixed the staggered moment m s by its Hartree Fock 
value, because for the doping values under consideration 
no long range order occurs. The dispersion in this anti- 
ferromagnetic background is a clear continuation of the 
evolution of the paramagnetic state with decreasing dop- 
ing. Obviously, a clear precursor of the antiferromagnetic 
state occurs already in the paramagnetic phase without 
broken symmetry. 
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vector k = (it, ir/8), which is shifted by (ir,ir) with re- 
spect to the main FS, the shadow band crosses the Fermi 
level. The intensity of the shadow states (ca. 10% of 
the main peak) agrees with the experimental observa- 
tion [Q. The possible occurrence of such shadow states 
for a system without long range antiferromagnetic order 
was originally proposed by Kampf and Schrieffer p3| . 
Using a phenomenological ansatz for the spin suscepti- 
bility, they argued that for sufficient long ranged anti- 
ferromagnetic correlations a coupling of states k and k' 
with k' — k — Q| < might lead to distinct shadow 
states. From their calculation one can estimate the mag- 
netic correlation length £ to be of the order of 20 lat- 
tice spacings, to obtain observable satellite peaks. This 
is much larger than the correlation length observed i n 
neutron scattering experiments which is only a few 
lattice spacings. 




FIG. 2. Quasiparticle dispersion in the neighborhood of 
the X-point for different doping concentrations. Note that 
for smaller doping (n = 0.44) two bands exist near the Fermi 
energy in the paramagnetic state. The dispersion of the main 
band (solid squares) resembles that of the highly doped sys- 
tem (dashed line). Considering both, the main band and the 
shadow band (open-crossed squares) one recognizes the evo- 
lution towards the two branches of the dispersion in an anti- 
ferromagnetic background (solid lines). 

In Fig. ||, we show the occurrence of shadows of the FS 
by plotting the spectral density £>k(w) for k-points be- 
tween the X and M-point. Besides the main peak, not 
crossing the FS, satellite peaks near the Fermi energy 
can clearly be observed. These satellite peaks are the 
new quasi particle states of Fig. (open-crossed squares), 
which build up the shadows of the FS. For the wave 



FIG. 3. Spectral density for k-points between the X and 
M— point, i.e. k = (tv, ^ir) with I = 4 • • • 10. The main peak 
reflects the dispersion of the bands similar to the case for 
U — 0. The satellites, indicated by the arrows, build up the 
shadows of the Fermi surface. The inset shows the imaginary 
part of the self energy for different momenta. 

Therefore, it was argued in Ref. Q] that there can- 
not be a well defined kinematics of the antiferromagnetic 
spin excitations that lead to the observations of Aebi 
et al. However, in our calculation the correlation 
length, obtained from the k-dependence of Vk(o;), is 2.5 
lattice spacings, in agreement with those experiments. 
Furthermore, we find that for larger doping (n = 0.4) 
the shadow states at low energies vanish, although the 
correlation length only decreases slightly. In order to 
resolve a shadow peak at k' = k + Q, the exitation ener- 
gies of the states near k should be within an energy range 
smaller than the inverse lifetime of the state k + Q, i.e. 
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I £ k— Ek+c- 1 ! < |ImSk+Q|- This can be fulfilled for a large 
£ or for a band with large effective mass. Because of the 
small correlation length of our calculation, the shadow 
states can only be understood from the above discussed 
strong deformation of the quasi particle dispersion. 

In the inset of Fig. 0, we show the frequency de- 
pendence of ImSk(^) for various k-points. We find a 
pronounced difference between momenta at the FS and 
away from it. At the FS, the transition from a low 
frequency io 2 to a linear in ui behavior occurs below 
0.005 eV. This extremely low energy scale is expected 
to be of importance for the anomalous transport prop- 
erties of the high-T c materials. At k = (tt, 7r/8), where 
the shadow band crosses the Fermi level, we find a dou- 
ble well structure. This reflects the strong coupling of 
the shadow states to the main FS and is a precursor of 
the singular behavior of ImSk(w) in the antiferromag- 
nctic state This anomalous frequency dependence 

is different from the main assumption of the Luttinger 
theorem 22 (ImEk(w) oc ui 2 ) and leads to an occupa- 



tion of shadow states and consquently to a violation of 
the theorem, shown in Fig. [lj. An experimental indication 
for the violation of the Luttinger theorem in underdoped 
systems was recently found by Liu et al. 0. 

In conclusion, we presented a new numerical method 
for the solution of the FLEX-equations on the real fre- 
quency axis, which permits the analysis of the fine struc- 
ture of the excitation spectrum for high-T c materials. 
A central role for all results shown here plays the ex- 
ceptional behavior of the quasi particles at the Appoint. 
The occurrence of a fiat dispersion at the X-point was 
shown to be of extreme importance for a deep under- 
standing of the shadows of the Fermi surface without 
broken antiferromagnetic order and for the violation of 
the Luttinger theorem. All this shows clearly that the 
quasiparticle properties of the high-T c superconductors 
are closely intertwined with the strong antiferromagnetic 
fluctuations. In view of the importance of antiferromag- 
netic correlations for superconductivity jjj], we extended 
our approach to the superconducting state, to investigate 
the phenomena discussed in this paper below T c . |16| 
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